Atomistic, microstructural and micromagnetic aspects of the multiscale modelling of 

hysteretic phenomena 
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We formulated a technique which combines the first-principles, micromagnetic and microstruc- 
tural calculations and allows us to study the nature of hysteretic phenomena in hard magnets. 
Two distinct sources of coercivity in polytwinned CoPt type magnets, domain wall pinning at an- 
tiphase boundaries and splitting at twin boundaries, are illustrated for a realistic microstructure. 
Methodology of multiscale modelling of hysteretic phenomena in nanoscale magnets is discussed. 

PACS numbers: 75.60.Ch, 73.20.-r, 75.50.Ww, 75.60.Jk 



I. INTRODUCTION 

It is a well established experimental fact that the hys- 
tetic properties of a magnetic material are very sensetive 
both to its magnetic atomic constituents and to its mi- 
crostructure. The studies of pure magnetic phenomena 
on the nanosize scale (domain wall structure, shape of 
domains and so on) can usually be addressed with appro- 
priate micromagnetic techniques (e.g., the finite element 
method). Quite extended regions (1-1000 nm and up) 
can be studied with these techniques, but such important 
ingredients of the problem as the atomic-scale properties 
of defects and the microstructure are usually addressed 
only phenomenologically. The first problem is related 
to the fact that no first-principle calculations have been 
used to provide adequate description of interatomic in- 
teractions on the atomistic scale, whereas the second one 
is due to lack of adequate theoretical description of the 
microstructure. This microstructure plays an especially 
important role in the properties of hard magnets, because 
their high coercivity always develops in quite characteris- 
tic nanoscale microstructures containing a high density of 
such microstructural elements as grain boundaries, twins, 
interphase and antiphase boundaries (APB), etc. Al- 
ready here it is evident that the theoretical description 
of the effect of microstructure on the hysteresis loop of a 
hard magnet presents a rather complex task due to the 
inherently multiscale nature of the problem and the pres- 
ence of several entirely different interactions. Interaction 
of domain walls (DW) with many important defects is 
determined by variations of the microscopic interaction 
parameters within the regions of atomistic size (~1 nm). 
The DW width (5-10 nm in most hard magnets) is an- 
other length scale, while the microstructure itself has one 
or more additional length scales (typically within the 10- 
200 nm range). Each of these length scales is physically 
important, and all of them must be linked together in 
order to describe the magnetic properties consistently. 
Therefore, a simultaneous inclusion of magnetic interac- 
tions on different scales (from 0.1 nm to 1000 nm) and 
corresponding elastic interactions (responsible for a given 
microstructure) seems to be important for the consistent 
description of the hysteretic phenomena. 

Whereas the qualitative considerations above have 



been around for a long time, to our knowledge, no the- 
oretical methods combining the atomistic, microstruc- 
tural and micromagnetic parts of the multiscale coerciv- 
ity problem in real systems have been suggested. Be- 
low we show how one can combine the first-principles, 
micromagnetic and microstructural calculations and ap- 
ply them to study the rich physics of hysteretic phenom- 
ena in hard magnets of CoPt type and demonstrate that 
coercivity has a natural multiscale character and must 
be studied using an approach including all the relevant 
scales and interactions on an equal footing. 

Below we first discuss the choice of studied CoPt type 
systems in section II. Microstructure of these materials 
and its theoretical description is described in section III. 
In section IV we outline the microscopic mean-field tech- 
nique replacing the continuous micromagnetic approach 
for the studies of microscopically defined microstructures. 
The results on atomistic scale obtained in ab initio calcu- 
lations will be discussed in section [v] in connection with 
APB defects. In section VI we combine the description 
of all the relevant length scales and interactions and find 
the structure of a macrodomain wall in a realistic poly- 
twinned microstructure. Section VII concludes the pa- 
per. 



II. CHOICE OF SYSTEMS 

Intermetallic hard magnets CoPt, FePt and FePd de- 
velop high coercivities in the tetragonal Llo phase. All 
the accumulated knowledge about this magnet family 
suggests that it can serve as a prototype for the prob- 
lem we specified in introduction. First of all, the main 
microstructural features in these alloys are well estab- 
lished experimentally on the nanoscale level (10-100 nm) 
which is comparable with S. Notably, the microstruc- 
ture is dominated- by coherent crystallographic defects 
(twins and APB'sEJ) and not by grain boundaries, as in 
other hard magnets. This circumstance provides two ma- 
jor advantages: the microstructure may be consistently 
simulated theoretically (see section III), while the mi- 
croscopic properties of defects are tractable in modern 
first-principles calculations. Indeed, the description of 
electronic structure of alloys of 3d atoms with Pd and 
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Pt is sufficiently reliable, and the range of perturbation 
near an APB (1-5 nm) may also be covered. Such pa- 
rameters as magnetic anisotropy and effective exchange 
coupling are also routinely computed. With this method- 
ical background, identification of the links between the 
microstructure and magnetic properties of bulk CoPt- 
type magnets is of interest both for the understanding 
of their own physics, and for the general theory of hys- 
teretic phenomena. Knowledge of such links may prove 
useful in the design of future nanofabricated magnets. 



III. MICROSTRUCTURE: POLYTWINNING 
AND ANTIPHASE BOUNDARIES 

Typical thermal processing for CoPt type alloys in- 
volves high-temperature annealing in the disordered Al 
(fee) phase area of the phase diagram followed by a rapid 
quench and ageing at a lower temperature in the Llo 
phase area. It is known from [transmission electron mi- 
croscopy (TEM) experiments!™ that the microstructural 
evolution during ageing includes two stages, the 'tweed' 
stage, and the twinning stage. The 'tweed contrast' cor- 
responds to a uniform pattern formed by small (1-10 nm) 
ordered domains with a dominant {101} orientation of in- 
terfaces and essentially random distribution of the c-axis 
directions. When the growing ordered domains reach 
some characteristic size (10-20 nm), a new pattern de- 
velops, with the formation of large 'polytwinned' stacks 
containing ordered 'bands' ('c-domains') with two alter- 
nating directions of the c-axis making tt/2 angles with 
each other. The interfaces between the c-domains within 
a polytwinned stack are all parallel to each other and 
lie in one of the 101 planes. The main features of mi- 
crostructure on both these stages stem from the presence 
of three possible orientations of the c-axis of the ordered 
phase relative to the cubic disordered phase. 

Physically, polytwinning is explained by the fact that 
such structures eliminate the volume-dependent part of 
the elastic energy stemming from the coexistence of or- 
dered domains with different c-axis directionstla Main 
features of the microstructural evolution at tweed and 
polytinning stages were reproduced using a phenomeno- 
logical 2D modelcl. Recently a microscopic model for the 
description of elastic interactions in alloy&i undergoing 
Llo type ordering was proposed and usedB in detailed 
studies of various aspects of microstructural evolution. 
In such alloys the elastic interaction turns out to be ef- 
fectively non-pairwise, but it retains the standard long- 
wavelength 'elastic singularity' of the k~ 2 typeu. There- 
fore, the elastic interaction, as the dipole-dipole inter- 
action, formally has an infinite range. However, due to 
the fact that the elastic energy of a non-polytwinned ar- 
ray of ordered domains grows faster with the domain size 
compared to the contribution from the interface tension, 
there is a characteristic size lo at which the elastic in- 
teraction begins to affect the microstructural evolutionH. 
This size (usually about 10 nm) corresponds to the aver- 



age size of domains on the tweed stage; the thickness of 
twins on the polytwinnecLstage can not be less than Iq. 

As it was shown earlieru, the twin boundaries strongly 
affect the structure of 'macrodomain walls' (MDW) 
crossing the polytwinned stacks by splitting them in seg- 
ments that are coupled only magnetostatically. Each seg- 
ment is pinned individually by pinning centers within the 
c-domain. This leads us to another generic feature of 
the polytwinned stage — the c-domains are interspersed 
with antiphase boundaries (APB) between ordered do- 
mains with the same direction of the c-axis and an an- 
tiphase shift along this axis. The APB's act as pinning 
centers for DWs due to a local suppression of MCA (see 
Section |v|). Clearly, the efficiency of pinning depends 
on the pattern formed by APB's. So far, very little is 
known about the dependence of this pattern on the al- 
loy composition and thermal processing. For example, 
in Fig. 5 of Ref. || the APB's in a CoPt sample appear 
to be rectilinear, while the available TEM imagesa for 
FePd alloys reveal quite dense patterns of curved APB's, 
albeit with some tendency to a preferred alignment. The 
former case is readily reproduced in microscopic simula- 
tionsQ, while the prerequisites for the latter pattern are 
unclear. It is quite obvious that pinning efficiency should 
strongly depend on the prevailing APB pattern, and each 
case should be studied separately. One way to under- 
stand the combined effect of MDW splitting and pin- 
ning in various microstructures is to take characteristic 
microstructures obtained in microscopic simulations and 
add the magnetic degrees of freedom using the method 
described in section IV. Since magnets operate at rela- 
tively low temperatures where diffusion is suppressed, the 
microstructure may be fixed in magnetic simulations. In 
this paper we will illustrate this approach by finding an 
equilibrium structure of MDWs in a polytwinned mi- 
crostructure taken from a microscopic simulation where 
APB's are predominantly rectilinear. 



IV. MAGNETIC SIMULATIONS: 
MICROSCOPIC MEAN-FIELD METHOD OR 
MICROMAGNETISM 

The simplest way to add the magnetic degrees of free- 
dom to the microscopically simulated microstructure is to 
use the atomistic version of the micromagnetic approach 
which turns out to be convenient in our problem where 
the relevant length scales are of the order of 5-50 nm. 

Consider a binary alloy AB with the classical Hamil- 
tonian 

H = H COI a{m} + UiUj f— Jijpifij + fl t D i3 flj 

i<j 

+ y / n i {e i (fli) - H /Zj] (1) 

i 

where H con { is the configurational part of the Hamilto- 
nian; i and j run over lattice sites; rij = 1 if site i is 
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occupied by a magnetic atom A and = otherwise 
(let B be non-magnetic for simplicity); jli is the rigid 
classical magnetic moment of the atom at site i; Jij, the 
exchange parameters; Ho, the external magnetic field; 
ei(pi), the MCA energy equal to — &j(7^ej) 2 for easy-axis 
anisotropy; and Dy, the magnetic dipole-dipole interac- 
tion tensor. 

Due to the fact that exchange interaction dominates 
at small distances, in the vast majority of problems the 
magnetic moments may be considered to be everywhere 
close to ferromagnetic alignment, and the angles between 
/Vs with non- vanishing Jy are usually relatively small. 
Therefore, the exchange term in H should be considered 
as the expansion of the total energy in the angles of devi- 
ation from this alignment. The parameters Jij are thus 
identified with the second derivatives of the total energy 
in respect to these angles, rather than with the parame- 
ters of the Heisenberg model defined to account for arbi- 
trary angles between fits. 

The inhomogeneous states of the system may be de- 
scribed by the free energy 

F=(H + T\nP) (2) 

where T is temperature, and P is the distribution func- 
tion. For the canonical statistical ensemble InPo = 
f3(F — H), and the system is homogeneous. For inho- 
mogeneous magnetic states with a fixed atomic configu- 
ration, let us introduce the "generalized Gibbs distribu- 
tion" P — exp[/3(F — H c g)] as it was donecl for H con {. 
In the theory of second-order phase transitions the effec- 
tive Hamiltonian iJ ft-|is a functional of the macroscopic 
order parameter fielcfl, while in our microscopic case it 
contains some 'effective parameters' depending on the 
environment, i.e. on the average values of jli and, if nec- 
essary, of their powers and correlators. In the simplest 
mean-field approximation (MFA) which will be used in 
this paper, in equilibrium these parameters are the 'mean 
fields' Hi = H -I- Ylj{Jij ~ Dij)cjUij where m 3 = 
and Cj — (rii) are the local magnetizations and concen- 
trations. Specifically, in this approximation we have 

H cS = 2J «» [-Hj/tj + ei(Jk)] , (3) 

i 

and the free energy 

F=(H-H cS )-T In Tr exp(-f3H cS ) (4) 
after averaging assumes its final form 

F = -Ej. DD -T^Cjhx J dpi exp [/3(H^ - e»)] . (5) 

i 

Here Ej^dd — J2 c i c j m i(~Jij + ^ij) m j ^ s the total ex- 
change and dipole-dipole energy. Equilibrium states for 
given boundary conditions may be found by solving the 
'self-consistency' relations mi = —dF/dHi. 

Let us emphasize that the above form of P is actually 
a reduced Gibbs distribution. The effective Hamiltonian 



contains effective parameters depending on the averages 
over dynamical variables, while the Gibbs distribution 
contains the Hamiltonian with actual interaction param- 
eters. In particular, the generalized Gibbs distribution 
with the effective Hamiltonian and proper boundary con- 
ditions may be used to describe inhomonegeous equilib- 
rium states (e.g., the DW's), while the standard Gibbs 
distribution may not, because it retains the full symme- 
try of the Hamiltonian. 

In cases when MFA is inadequate (e.g. in frustrated 
systems) one can use more refined statistical methods to 
calculate F, e.g. the-cluster variation methodEj or the 
cluster field methodcil, although this will considerably 
complicate the calculations. 

The microstructure is defined by the set of Cj. For the 
studies of simple configurations this set may be prepared 
'by hand'. More realistic configurations may be obtained 
in microscopic simulations!! 

If magnetization M(r) slowly varies in space and is 
constant in magnitude, Eq. (^) reduces to the micromag- 
netic free energyEH In this case all choices of Jij and bi in 
the defect-free regions are equivalent if they produce the 
same macroscopic properties A and K. However, varia- 
tion of Jij and bi at interatomic distances near defects like 
APB's must be studied using first-principles techniques. 
MFA (or more refined) calculations with these parame- 
ters may be used to describe DW interaction with the 
defect at the length scale of S. At larger., microstructural 
length scales micromagnetic methoda3E-2 may be used 
with singularities of A and K at the defects. However, 
as we will see below, in hard magnets the microscopic 
approach also turns out to be convenient for the studies 
of regions containing up to ~ 10 6 atoms; in such cal- 
culations some model Jy and hi reproducing the actual 
defect properties may be used. 



V. ATOMISTIC SCALE: PINNING AT 
ANTIPHASE BOUNDARIES 

As it was noted above, the DFT methods are quite 
reliable for d-systems like Co-Pt. These methods were 
not used to study the effect of lattice imperfections on 
the hysteretic phenomena because the extended range 
required to simulate a DW is certaily inaccessible for any 
band structure technique. However, for CoPt-type alloys 
such study is feasible both in a direct ab initio approach 
(large supercell modelling of DW and defect) and on the 
model level (calculation of parameters of Hamiltonian). 
In the present paper we are setting up a technique com- 
patible with microstructural simulations, so we will dis- 
cuss only the latter approach. Before we proceed with 
this calculation let us first discuss the physics of hys- 
teretic phenomena on the atomistic level. 

As we noted above, polytwinned alloys usually contain 
a high density of APB's in the c-domains; it was ar,«pjl 
that pinning on them may explain high coercivityl 15 !! 16 !. 
In this section we will show using first-principles calcu- 
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lations that pinning at APB's (developing on atomistic 
scale alone) is indeed an important source of coercivity 
in given alloys. 

If pinning is generated by planar faults like APB's, 
the maximum unpinning threshold H u is achieved when 
DW segments (DWS) are parallel to the faults. If the 
planar fault is represented micromagnetically as a slab of 
thickness w <C 8 with modified exchange and anisotropy 
constants A' and K', then H u is given byc3: 



the free energies of the two large samples without and 
with the defect. Thus, we obtain 



~H~a 

3/2 ^ 



A 



K>_ 
~K 



(6) 



0.60 (our definition of 5 contains 



with a — 7r/3 
an extra factor tt compared to the notation of Ref. 17). 
Positive or negative H u corresponds to DW attraction to 
or repulsion from the fault. 

One should bear in mind that formula (j^) may be 
used only if the range of exchange interaction (distance 
at which A converges) is smaller than 8; otherwise, the 
representation of the planar fault by any variation of A 
across the fault is meaningless. In fact, this is the case in 
all hard magnets, because A is proportional to the sum 
J2j Jij( r ij n ) 2 ( n is the unit vector normal to the APB) 
which converges very slowhU. Contrary to the homoge- 
neous case, A(x) can not be conveniently computed using 
a summation in the reciprocal space (x = rn is the co- 
ordinate normal to the APB). Therefore, one should find 
the 'exchange contribution' to H u directly using the mi- 
croscopic representation of the defect and the formalism 
of Section IV. Note that the problem of slow conver- 
gence occurs only in the continuous formulation, because 
the above second-moment sum applies only to the re- 
gion where the DW profile m(x) may be replaced by its 
second-order Taylor expansion. The parameters J« may 
be calculated using the local spin density functional ap- 
proach^ which gives the second derivatives of the total 
energy over the deviations of a pair of spins from the 
equilibrium ferromagnetic alignment, in accordance with 
our definition of J^ . Technically, H u may be found as 
the maximum field at which an equilibrium state of the 
DW at the APB is possible. 

The modification of K at the defect has a negligi- 
ble effect on the DW profile m(ir) and hence may be 
accounted for perturbatively. Therefore, we may di- 
vide H u in two parts, H u = H ue + H ua , where H ua — 
aH a (w/S)(l — K'/K), and H ue must be found micro- 
scopically as described above, assuming that MCA is not 
affected by the defect. 

It is convenient to rewrite H ua in a different form more 
suitable for first-principles calculation. Due to the fact 
that the MCA modification at the defect may be treated 
perturbatively, the above formula for H ua does not de- 
pend on how MCA behaves close to the defect, as long 
as the perturbation falls off at distances smaller than 
8. In the 'slab' model used to derive (0) the product 
A/ a = w(K — K') may be identified as the 'anisotropy 
deficit' per unit area of the defect, i.e. the difference in 



H ua = aH a Af a /K5. 



(7) 



The value Ae a of A/ a at T = may be directly found 
using the first-principles calculations. To this end, it is 
necessary to take a sufficiently large supercell containing 
the defect, calculate its MCA energy, and subtract the 
corresponding MCA energy of the pure material. In this 
context, 'sufficiently large' means that Ae a is converged 
with the supercell size. One should bear in mind that 
symmetry is lowered at the defect, and MCA modifica- 
tion at the defect may lack the easy-axis symmetry of the 
pure material. 

We studied the modification of exchange and MCA 
at an isolated (lOl)-oriented APB in CoPt, FePt and 
FePd using the tight-binding linear-muffin-tin-orbital 
(TB-LMTO) method. We used elongated supercells of 
up to 28 atoms with an APB in the middle (the geom- 
etry of the (101) APB allows one to introduce a single 
APB per unit cell of 8n + 4 atoms). 

The analysis of the exchange parameters for CoPt 
in such supercells in conjunction with MFA simulations 
shows that changes in the local geometry together with 
the modification of Jy at such APB generates only a 
weak DW repulsion from the APB with \H ue \ < 1 kOc. 
The total exchange interaction between the pairs of 
atomic layers parallel to the (101) APB is notably af- 
fected only in the close vicinity of the APB, particularly 
for the next-nearest layers across the APB where extra 
nearest neighbors appear compared to the pure material. 

From the other hand, we found that MCA is strongly 
suppressed at the (101) APB leading to a considerable 
DW attraction to it. The values of Ae a per area of the 
APB corresponding to one formula unit in each atomic 
layer are 5.0E a for CoPt, 2.8E a for FePt and 4.8£ a for 
FePd, where E a is the bulk MCA energy per formula 
unit. Since high MCA is associated with Llo ordering, it 
is natural that local disorder introduced by an APB sup- 
presses MCA; one may expect it to be a generic feature 
of APB's and other defects. 

At low temperatures compared to the ordering transi- 
tion the order parameters and atomic structure of defects 
are almost independent on temperature. In this tempera- 
ture range (going up and beyond the room temperature) 
it is reasonable to expect that A/ a and K depend on T 
in a similar way. Therefore, an accurate estimate of H ua 
may be obtained using (^) with the zero-temperature val- 
ues of Af a and K and actual values of H a and 5 at the 
given T (the dependence of 8 on T is also quite weak, 
because both A and K decrease with increasing T). At 
higher temperatures the disorder at APB's may increase 
considerably, leading to a different T dependence of A/ a . 

Using (m) and the values of H a and 8 from Ref. 20 
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obtain the threshold fields of 6.8, 3.7 and 1.2 kOe for 
CoPt, FePt and FePd, respectively; these values are con- 
sistent with observed coercivities. Disorder at an APB 
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(which is always present in real allojjs with non-nearest- 
neighbour configurational interactions) notably increases 
Ae a , and hence H ua . This effect should also be expected 
at high temperatures affecting the temperature depen- 
dence of the coercive force. 



VI. MACRODOMAIN WALLS IN A REALISTIC 
MICROSTRUCTURE 

Although experimental observations are too scarce for 
reliable conclusions, in actual microstructures the APB's 
are rarely rectilinear, except for advanced stages of an- 
nealing (like the sample shown in Fig. 5 of Ref. ||); at 
early stages of annealing the c-domains contain a large 
density of intricately curved APB'sa. Therefore, H u in 
real alloys may be significantly lower compared to the val- 
ues obtained above for the plane-parallel configuration of 
DW's and APB's. The best coercive properties should be 
achieved as a result of the competition of a number of fac- 
tors depending on the annealing time, most notably H u , 
H s , Id and l m . A complicated interplay of these factors 
may be studied by direct simulation of MR in realistic 
microstructures obtained in microscopic simulations. To 
illustrate this approach, we found a stable configuration 
of two MDW's in the CoPt model described above with a 
microstructure similar to the experimental one shown in 
Fig. 5 of Ref. 0. We took a simulated structure similar to 
Fig. 5 of Ref.Tl but obtainedEil with a larger tetragonal 
distortion e m — 0.2, and augmented it to a 256x256x1 
bet simulation box (~ 70 nm) using the periodic bound- 
ary conditions. In order to reproduce DW attraction to 
APB's in accordance with our band structure calcula- 
tions, we put bi = in the close vicinity of APB's (de- 
fined by a certain combination of local order parameters), 
and hi — b at all other sites with easy axis coinciding 
with the local tetragonal axis (also found using local or- 
der parameters). The exchange parameters were set 
as it was specified above, making the exchange constant 
A continuous at (101) APB's. This choice of bi and Jy 
follows the general prescription given above after Eq. (||) ; 
the unpinning threshold for parallel DW and APB was 
made to be close to its value obtained in Section |v| for 
CoPt. 

The resulting magnetic structure is shown in Fig. [I] il- 
lustrating and extending our conclusions about the prop- 
erties of MDW. The MDW's are seen to be 'soft'; DWS's 
in the c-domains have a considerable freedom in their 
relative displacement, and several of them are pinned at 
APB's. Notably, the DWS's do not change their shape 
strongly when partially pinned by APB's (e.g., in the 
upper-left corner). Strong magnetic fields 'left behind' 
by displaced DWS's are responsible for the attraction of 
DWS's in an isolated MDW. 

Two separate sources of coercivity in polytwinned 
CoPt type magnets, domain wall pinningp-at antiphase 
boundaries and splitting at twin boundaries^, are clearly 
present in our simulations for a realistic microstructure. 
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FIG. 1: Macrodomain walls simulated for the model of a poly- 
twinned CoPt magnet (quasi-2D box with ~ 70 nm edge). 
The distributions of magnetization (top) and dipole fields 
(bottom) are virtually planar and are shown by arrows; color 
shows the magnetic charge density p = — divM (red for posi- 
tive, blue for negative). The black lines show horizontal twin 
boundaries and APB's within the c-domains. 



Splitting may strongly affect the efficiency of DW pin- 
ning, and it is important to account for both these effects 
simultaneously in order to describe the coercivity. 



VII. CONCLUSION 

We have described the DW properties in polytwinned 
CoPt-type systems where two sources of coercivity, do- 
main wall pinning at antiphase boundaries and splitting 
at twin boundaries, are included simultaneously using a 
realistic microstructure. Pinning is related to MCA low- 
ering at the planar defect (APB) whereas splitting comes 
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from the competition between the long-range dipole- 
dipole and short-range exchange interactions. 

From the applied point of view, it is usually as- 
sumed that polytwinning is detrimental for high coer- 
civity. However, splitting may increase the efficiency of 
DW pinning by allowing the DWS's to 'find' pinning sites 
independently. It seems that the application potential 
of polytwinned microstructures has not been exhausted, 
and better properties may be achieved on the basis of 
this new theoretical understanding. 

From the methodical point of view, we have shown 
that the coercivity is an intrinsic multiscale property in 
hard magnets and to describe this property one has to 
use a theory which necessary combines a synergy of dif- 
ferent length scales in a presence of physical interactions 
of entirely different nature. Whereas the different aspect 



of the problem of magnetization reversal can be stud- 
ied using one level approach, the problem as a whole 
requires more complicated description. We believe that 
such synergistic type of approach represent a generic fea- 
ture for any theory to be succesful in area of nanoscale 
type of magnets: bulk systems, multilayers, clusters and 
nanowires. 
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